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Abstract 

New cross sections for the reaction ep e'rjp are reported for total center of mass energy W=1.5- 
2.3 GeV and invariant squared momentum transfer Q^=0. 13-3.3 GeV^. This large kinematic range 
allows extraction of new information about response functions, photocouplings, and rjN coupling 
strengths of baryon resonances. A sharp structure is seen at W ^ 1.7 GeV. The shape of the 
differential cross section is indicative of the presence of a P-wave resonance that persists to high 
Q^. Improved values are derived for the photon coupling amplitude for the 511(1535) resonance. 
The new data greatly expands the range covered and an interpretation of all data with a 
consistent parameterization is provided. 
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I. INTRODUCTION 



Photoproduction and electroproduction experiments on the nucleon provide a clean probe 
of nucleon structure since Quantum Electrodynamics is well understood. As a result, the 
matrix elements for N*,A* transitions, commonly called the photon coupling am- 

plitudes, are sensitive to the nucleon and A^* quark-level wave function. These amplitudes 
have traditionally been calculated using quark models l[ 2^ Q], but recently progress has 
been made in applying the techniques of Lattice QCD 6|. Experimental measurements 
are currently being made of a number of different baryon resonances in several different final 
states. For a review of the current status, see Reference 0, Q]- 

Disentangling the wide and overlapping states that populate reaction data has been a 
long-lasting problem. In the mass region above total center of mass (cm.) energy (W) of 
1.5 GeV, many overlapping baryon states are present and some are not well known. The 
reaction ep — > e'rjp is especially clean, since processes involving rjN final states couple only 
to isospin | resonances, simplifying the analysis. A prominent peak in the total cross section 
is seen for rj production at W = 1.535 GeV in both 7A^ and vrA^ experiments. This is widely 
interpreted as the excitation of a single resonance, the spin |, negative parity, isospin | state 
5*11 (1535) ^ This state has a branching ratio to rjN of 45-60% compared to at most a few 
percent j9|, [lOj for other states. This is a very interesting and unusual pattern. 

Eta photoproduction experiments have reaffirmed the strong energy dependence and S- 
wave (isotropic) character close to threshold 11|. Using polarized photons [3], new values 
for T]N decay branching ratios of other resonances have been determined through interference 
with the dominant S'ii(1535). 

Electroproduction cross sections can be used to extract the photocoupling amplitude for 
non-zero values of the squared momentum transfer (Q^) from the electron to the resonance. 



Using rj electroproduction [13I, 



14 



16 



17 



18 



19j, an unusually flat dependence of 



the photocoupling amplitude was found for the 5*11(1535) in contrast to the nucleon form 
factors and photon coupling amplitudes of other established resonances, e.g. P33(1232) [20|. 
Although previous 1] angular distributions were largely isotropic at all Q^, no detailed re- 
sponse functions were extracted because of the poor angular coverage using traditional 
magnetic spectrometers. Although the dependence was clearly different than for other 
resonances, the results were comprised of many different experiments whose results appeared 
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to be inconsistent with each other. An analysis by Armstrong et al. [19] showed that much 
of the inconsistency was due to different assumptions about 5*11(1535) properties used by 
the individual experiments. 

In our previous publication {21] we presented results on rj electroproduction based on the 
first data taken with the CEBAF Large Acceptance Spectrometer (CLAS) 22j at Jefferson 
Lab. We extracted the photocoupling amplitude (^1) for the 5*11(1535) over the range 0.25 
GeV^ < <1.5 GeV^ from our data. In addition, we observed indication of a structure at 
W ^ 1.7 GeV in the total cross section which is also seen as a change in the shape of the 
differential cross section at the same energy. 

The energy region around W ^ 1.7 GeV has received significant attention lately. A 



CLAS TT+TT electroproduction experiment 231] found excess strength at the same energy 



beyond theoretical predictions based on previous data. This excess strength was tentatively 
identified as a P-wave resonance; either the decay properties of Pi3(1720) change significantly 
or there is a new spin |^ state. A recent r] photoproduction experiment at Bonn 24] provides 
a comprehensive set of cross section data from near threshold to well beyond the resonance 
region. In the same paper 2J] a partial wave analysis of these and other r] electroproduction 
data finds strong excitation of a J = 3/2~^ state at 1775 MeV which they identify with 
the Pi3(1720). Although previous analyses jo] found weak evidence for any rjN decay of 
resonances at W ~1.7 GeV, the new data is of much higher quality than the older data. 

The data and analysis reported here use a data set taken with the same apparatus (CLAS) 
as used in our first publication [21]. The new data have an order of magnitude more 77 events 
than in that previous paper and a much larger kinematic range. Therefore, the new values 
presented here supersede the previously published data. Our reach in (0.13 GeV^-3.3 
GeV^) is more than twice as large as in our first publication. This allows a large extension 

of the range where we can extract the photocoupling amplitude Ai of the proton to 

2 

5*11(1535) transition. We also more precisely determine the non-isotropies in the differential 
cross section and show evidence for a significant contribution to rj electroproduction due to 
a P-wave resonance with a mass around 1.7 GeV. 

The paper first presents some formalism needed to understand the measurement and its 
analysis, followed by details of the experiment. We then present results for the inclusive and 
exclusive analyses. Discussion of these results with a Breit-Wigner model and conclusions 
complete the paper. 
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II. FORMALISM 



The kinematics for the ep e'rjp reaction are shown in Figure [H It can be characterized 
in terms of the squared 4- momentum transfer between the electron and proton (— Q^) carried 
by the virtual photon (7t,), the invariant mass of the 7i, — p system (W), and the scattering 
angles of the final state rj in the rest frame of the 7?, — p system {6*, 0*). These angles 
are also the decay angles of the resonance in its rest frame. We use the superscript * for 
quantities evaluated in this frame. The five-fold unpolarized differential cross section for 




Scattering plane 



Reaction plane 



FIG. 1: The reaction is depicted in the center of mass system, where the resonance is at rest. The 
meson decay polar angle is defined relative to the virtual photon momentum, and the azimuthal 
angle is defined relative to the the electron scattering plane. 



the ep —>■ e'rjp process at a specific energy {E) may be expressed as the product of the 
transverse virtual photon flux (F^) in the Hand convention [251] and the cm. cross section 
for virtual-photoproduction of the prj pair: 



rf5 



cr 



dWdQ'^dn* 

The cross section for the virtual reaction 'y^p 
display the dependence on (p*. 



dn* 



(1) 



rip is written by convention to explicitly 



dfl: 



-{ivP — ^ tip) 



ctt + eciL + 

A/2e(l + e)(JLTCOS < 
eaTTCOs 20*, 



+ 



(2) 
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where e is the longitudinal degree of polarization of the virtual photon and is given by 

(3) 



2 /,T,xn-l 



l + 2(l + lj)tan'(| 

where q is the magnitude of the 3-momentum of the virtual photon and \1/ is the electron 
scattering angle. Since e is invariant under coUinear transformations, q and \1/ may be 
expressed either in the lab or cm. frame. The cross sections for transverse and longitudinal 
photon are represented by ax and a^, respectively. In addition, glt is a contribution 
due to interference between transverse and longitudinal amplitudes and cttt describes the 
interference between amplitudes for the two different transverse polarizations, either aligned 
or anti-aligned with the spin of the target proton. All four of these terms depend on W , Q^, 
and cos^*. 

In order to identify individual baryon resonances, the cross section should be decomposed 
into partial wave amplitudes. These amplitudes are most often labeled by the electromag- 



netic multipole notation [26|]. Multipoles are commonly labeled Ei±^ Mi±, and Si±, where 
I is the orbital angular momentum of the final rjp system and ± denotes whether the total 
angular momentum is / ± |. E and M refer to electric and magnetic transitions involv- 
ing transverse virtual photons, while the longitudinal (S) transitions involve longitudinal 
photons. 

The response functions and multipoles have contributions from underlying resonant and 
nonresonant reaction mechanisms. When evaluated at the peak of the resonance, the mul- 
tipole is expressed in terms of both the photocoupling amplitude and the hadronic decay 
properties of the individual resonances in a commonly accepted way [^^j . The photocoupling 
amplitudes are labeled by the total helicity (| or |) and the virtual photon polarization 
(transverse or longitudinal) and depend on the invariant squared momentum transfer to the 
resonance (Q^). The shape of the resonance determines the W dependence of the resonant 
part of the multipole. Spin-^ resonances will be described by one transverse amplitude 
{Ai) and one longitudinal amplitude (5*1). In terms of multipoles, an 5*11 resonance has 
an Eq^ (electric dipole) and a 5*0+ transition; a Pu has Mi_ and Si^ transitions. Higher 
spin resonances will be described by both Ai and Ai photocouplings (and both E and M 
multipoles). Extraction of the multipole amplitudes from the cross section data, see e.g. 



281], is not unique because more than one bilinear combination of multipoles have identical 



angular distributions. We therefore choose simplified methods (discussed below) to analyze 



the data. 

The differential cross sections can be calculated with a model of resonance produc- 
tion/decay and the nonresonant processes. This cannot yet be done from a fundamental 
field theory such as Quantum Chromodynamics (QCD). Instead, models are used that have 
parameters determined from data. The r^-MAID 29| model uses an isobar model 30|] to 
construct the cross section for rj photo- and electroproduction; parameters are fit by com- 



parisons with previous results ll|, |2l|, |3l|] . We have calculations from the MAID code for 



our kinematics. To further understand our data, we also do Legendre polynomial fits to the 
angular distributions. Both these results are described in Sect. IIVBI 

To analyze our angle-integrated cross sections, we make a further simplification which is 
possible because the 5*11(1535) resonance is dominant near threshold. Therefore, we ignore 
the nonresonant amplitude. If one can isolate the contribution of a single Su resonance to 
the Eq+ multipole, the cross section takes a simple form, 

da p;w 



nipK 



(4) 



where K = {W"^ — m^/{2mp) is the equivalent real photon energy, p* is the momentum of 
the outgoing rj in the 6*11 rest frame, and nip is the proton mass. The longitudinal multipole, 
5*0+, in principal contributes and we do not have the data to make the separation. However 



Sqj^ has been found to be small 



14 



16. 



17 



18 



09 



171] , and it was therefore ignored in previous analyses 
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2l| . An isobar model analysis of rjp, ir'^p, and vr+ri CLAS electroproduction 



data (4l| confirms the assumption of a small longitudinal component. In this analysis, the 

value of Si/Ai is about 15-20%: this translates to a few percent contribution to the cross 

2 2 

sections measured here. With the assumption that a single resonance dominates the cross 
section and 5*1 is small, Ai for 7p —>■ 5*11(1535) can be determined from 27] 



(5) 



where Eq+ refers only to the contribution from the resonance which is evaluated at the peak 
of the resonance. If there are other contributions to -Eo+, a model is needed to extract 
the resonance contribution. This formula contains terms related to the final state decay 
of the S*ii: (the total width of the S*ii(1535)), 6^ (the branching fraction into the rjp 
final state), and p* are all calculated at the mass of the 5*ii (Wr). Our current lack of 
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knowledge of these parameters leads to a model dependence in the extracted values for Ai. 



This prompted Benmerrouche et al. 
less model dependence 



32l | to propose using a quantity for each resonance with 



mlKhr, 



^1 = W "I " A.. (6) 

A\ depends on the matrix element for the initial state N* transition, while 

is proportional to the product of the matrix elements for the — > A^* and A^* rjp 
transition. For 5*11(1535), = ^/2TTlm{Eo+{WR)). Although is more closely related 
to experimental values, Ai is more easily determined from calculations, e.g. using quark 
models. Whichever quantity is used, the model dependence still exists when comparing 
calculations to experiment. 

We use the same resonance parameterization in all our calculations. The relativistic 



Breit-Wigner form is taken from previous rj photoproduction work 33|, |3J] and extended to 
non-zero angular momentum as: 

^W = p|r|„,(»')r P) 

where a is a constant that contains the photon coupling amplitude and kinematic factors, q* 
is the photon 3-momentum in the resonance rest frame, is the partial width for A^* rjp 
decay, and F^t is the total width, 

Tfi is the bare width and Bg{p*) is a Blatt-Weisskopf penetration factor 35(]. If i=0, this 



factor is equal to 1. The momentum ratios 



34| approximately account for proper phase space 



effects for the various final states (vrA^, rjN, and tcttN where the phase space factors for the 
ttttN final state are ignored). They are weighted according to estimates of the branching 
fractions to each final state. This form has been successful in matching data, but is not 
unique. 
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III. DETECTOR AND ANALYSIS 

The CLAS facility [22^ was designed for efficient detection of multi-particle final states. 
The data used for this measurement was taken in 1999 at electron beam energies of 1.5, 
2.5 and 4.0 GeV. A cylinder of liquid hydrogen was used as the target. Two different 
targets were used, 5.0 and 3.8 cm long. Toroidal magnet coils separate CLAS into 6 largely 
identical sectors, each covering roughly 54° in azimuthal angle (j) (with smaller coverage at 
smaller polar angle). Tracking drift chambers (DC) in CLAS measure angles and momenta 
of charged particles for lab polar angles in the range 8° < ^ < 142°. Outside the DC, 
scintillation counters (SC) provide time-of-flight measurements with which we can separate 
the charged hadrons into pions, kaons, and protons. For lab angles 9 < 48°, threshold 
Cerenkov counters (CC) and Electromagnetic Calorimeters (EC) distinguish electrons from 
charged hadrons. 

For this analysis, events were selected with an identified electron and proton. Since the 
momentum 4- vectors of the beam and target are known, the 4- vector for the putative rj can 
be determined from these two final state particles. A fiducial cut on these particles was 
applied to avoid the regions near the magnetic coils and the edges of the CC where the 
acceptance is changing rapidly. The momentum of the electron was required to be above 
400 MeV in order to be well above the trigger threshold. 

Cross sections were calculated as a function of and W for the angle-integrated data 
analysis and as a function of Q^, W, cos 6'*, and 0* for the differential data analysis. Cross 
sections are determined in a standard way by determining the yield in each of many bins, 
correcting for detector acceptance, and normalizing by the beam intensity measured with a 
Faraday cup and the calculated target thickness. 

As discussed in Sect. [Tll the extraction of resonance properties comes from an analysis 
of the cos 6'*, 0* distributions at specific values of W and Q^. Distributions of these vari- 
ables covered by the apparatus are determined by geometry. The large acceptance of CLAS 
guarantees almost complete coverage in cos 6'* and 0*. The beam energies of the experi- 
ment coupled with CLAS provided data at a wide range of Q'^ and W. We compare the 
kinematic range for the new experiment with that available for the previously published rj 
electroproduction data in Table [H 

Events were divided into separate kinematic bins as detailed in Sec. IIVI For each bin 
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Experiments 


W (GeV) 


Q2 (GeV2) 


cos 


e* 


0* (deg) 


Daresbury[13l 


1.51 1.55 


0.15 ^ 1.5 


not g 


iven 


not given 


Bonn [141 


1.51 ^ 1.56 


0.2 0.4 


-0.766 - 


0.939 


~0 


DESY[15l 


1.5 ^ 1.7 


0.22 1.0 


-1 - 


1 


^ 180 


DBS Y [161 


1.49 ^ 1.58 


0.6,1-0 


-1 - 


1 


15 ^ 90 


Bonn [IT] 


1.44 1.64 


0.4 


-0.643 - 


0.866 


-40 40 


DESY^ 


1.49 1.8 


2.0,3.0 


-1 - 


1 


^ 120 


JLab[19] 


1.48 ^ 1.62 


2.4,3.6 


-1 - 


1 


^ 360 


JLab[21j 


1.5 ^ 1.86 


0.375 ^ 1.5 


-1 - 


1 


^ 360 


this experiment 


1.5 ^ 2.3 


0.13 ^ 3.3 


-1 - 


1 


^ 360 



TABLE I: Summary of the kinematic ranges of previously pubhshed data compared to this experi- 
ment. This experiment is an extension of Ref. pl|]. Values given are the maximum ranges for each 
experiment. 



the T] yield was determined by fitting the distribution of missing mass recoiling against the 
outgoing e-p system. An example fit in one bin is shown in Figure O 

The fit is the sum of a signal at the 1] mass and a background function. The signal shape 
has a radiative tail and is corrected for experimental resolution; the background function 
is a polynomial. We use the data to determine both shapes. Both functions must then be 
modified by the geometric acceptance for this reaction because it has a rapid variation with 
respect to the kinematic parameters. This method is an extension of what was used in the 
previous CLAS data analysis [21]]. 

The shape of the signal was modeled in 2 steps to reproduce all features seen in the data. 
It is first described by a delta function at the 1] mass {m^) plus an exponential above 
representing the radiative tail, 

S{m) = (1 - f)5{m - rrir,) + fe{m - m^)e-"^"'^\ 

The fraction of events in the radiative tail (/), and a parameter describing the slope in 
the exponential (a) were determined for each W-Q^ bin from Monte Carlo generator events 



containing radiative effects 



36|, |37j. This signal shape was then convoluted with a Gaussian 



representing the experimental resolution to obtain the final signal shape (an analytic func- 
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FIG. 2: Sample missing mass {Mx) spectrum for ep epX. The bin shown is for 1^=1.535 GeV 
and Q^= 0.6 GeV^. The dashed hne (right scale) shows the acceptance that is calculated with a 
Monte Carlo program. The sum of the r] signal shape and the raw background function modified 
by the acceptance function is then fit to the missing mass spectrum for each bin. In the figure, the 
dot-dashed curve is the raw background function, D^^g from this fit, while the dotted curve shows 
those same values when multiplied by the acceptance function. The solid curve shows the full fit. 

tion) used to fit for the rj yield. In the fit to obtain the yields, only the magnitude of the 
signal and background functions were free. All other parameters were determined separately. 
The rms resolution for the missing mass peak ranged from 4 MeV in the low bins up to 
12 MeV at the highest Q^. The experimental 1] mass, experimental resolution width, /, and 
a were first fit to simple functions of and W for both data and Monte Carlo to smooth 
out statistical variations. The experimental r] mass was found to be within 1 MeV of the 
accepted value. Estimated contributions to the systematic uncertainty by these choices of 
parameters were evaluated in a later step. 

The background comes from tttt production. Although it has a smooth dependence on 
Mx, no models are available. We fit the background with a simple polynomial {Dhkg) in the 
missing mass {Mx) constrained to be zero at the kinematic limit {rrijnax) as required by the 
decreasing phase-space. 



where Am = rrimax ~ Mx and bo is the overall strength. One example of this function is 




(11) 
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shown in Fig. [2l At the highest beam energy, a shghtly more comphcated function was used. 
Both forms contained one parameter (Am') that was determined from our data by fitting 
to a polynomial in W. As with the peak shape function parameters, these fit parameters 
were included in the systematic uncertainty determination. 

The main structure in the background fit function comes from the variation of the ge- 
ometric acceptance of the detector as a function of missing mass. We found that proper 
modeling of this acceptance was very important. We determined this acceptance using a 
separate Monte Carlo program that generated ep — epX events with the X mass thrown 
randomly across the fit region. After requiring the scattered electron and proton to be in 
the fiducial volume of CLAS, we compared the generated and accepted events in order to 
calculate the background acceptance function. We multiply our simple background func- 
tion with the calculated acceptance function to obtain the final background used in the fit. 
Examples of all three of these curves are shown in Fig. [2j 

In a small number of bins where the cross section is low, statistical fluctuations in the 
background can lead to a best fit value for the number of r/'s that is negative. In this case 
we follow the suggestion of the Particle Data Group (PDG) ^ and report a negative value 
with error bars for the cross section. This provides sufficient information for constraints 
from these bins to be combined with nearby bins in comparing to theoretical predictions. 

Acceptance for the ep —>■ erjp reaction was calculated using a GEANT-based Monte Carlo 
simulation The event generator included radiative effects using the peaking approxi- 
33], and the cross sections have been corrected for radiation. When making a 



mation 
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major improvement in published cross sections, development of an appropriate event gener- 
ator is important. We use the data as a guide; the final cross section is dominated by S and 
P waves with the S'ii(1535) the dominant structure seen. An iterative procedure matching 
analyzed Monte Carlo to real data was used to develop the event generator. The same fitting 
procedure used on data was applied to Monte Carlo events to calculate acceptance. The 
acceptance has significant variation across the bins with a maximum value of about 60%. 
When approaching the kinematic limit, the acceptance falls off rapidly. At the higher values 
of W, the proton goes forward where there is a hole in the CLAS acceptance. This causes 
problems for 0^ ~ 180°. We only report results in bins where the acceptance is greater than 
3% and where it is not changing rapidly. 

A detailed study of potential sources of systematic uncertainty was made. Since the rj 
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peak shape and the background shape included various parameters, all were studied. The 
parameters were varied within the error bars determined in the fit for each cross section value. 
Additional tests were made for variations in particle identification and in the momentum 
scale. The momentum uncertainty arose from uncertainty in the details of both the magnetic 
field map and the alignments of the various tracking chambers. Sensitivity to momentum 
determination was largest close to threshold and falls off with increasing W. Since most 
of the rj events are produced near threshold, this source dominates the average systematic 
uncertainty. Cross sections were recalculated with slightly tighter fiducial cuts and this 
variation was considered as a systematic uncertainty estimate. A variety of underlying 
physics models were used for evaluating the systematic on the radiative correction: using a 
single S-wave resonance or two, varying the mass and width of a single S-wave, or including 
a P-wave resonance. The quoted systematic uncertainty on the radiative correction includes 
these effects, but is dominated by Monte Carlo statistics in the calculation. The total 
systematic uncertainty for each bin in W, Q^, and cm. scattering angles was the sum 
of all the components added in quadrature. The average total systematic uncertainty for 
the angle-integrated cross sections was 3.3%, 3.9%, and 7.1% for data at 1.5, 2.5 and 4.0 
GeV, respectively. The corresponding average estimated systematic uncertainties for the 
differential cross sections were 5.1%, 5.2%, and 7.6%. The breakdown by source for the 2.5 
GeV data (the set from which the largest number of data points come) is given in Table HTl 
The estimated systematic uncertainties for individual data points were seldom larger than 
the estimated statistical uncertainty. 



IV. RESULTS 



A. Angle-integrated Cross Sections 

To get the angle-integrated cross sections, the events were binned in W and Q^, as shown 
in Table UTIl The 1.5 GeV beam energy data covers the range from 0.13 to 0.4 GeV^, 
while the upper two beam energies cover from 0.6 to 3.3 GeV^. Each bin is labeled by 
its centroid. Results are tabulated in the CLAS database |39|. These cross sections are 
presented in Fig. [3l The prominent peak at W ~ 1.5 GeV is primarily populated through 
intermediate excitation of the S'ii(1535) resonance. Fits to a Breit-Wigner relativistic form 
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Sys Err Source 


Angle-integrated (int) 


Differential (diff) 


7 


0.03% 


0.06% 




0.4% 


0.7% 


f 


1.3% 


1.4% 


a 


0.1% 


0.1% 


DhkaiMx) 


0.1% 


0.1% 


fiducial cut 


0.6% 


2.3% 


radiative corr. 


1.0% 


1.0% 


momentum scale 


2.6% 


4.4% 


total 


3.9% 


5.2% 



TABLE II: Summary of systematic uncertainties for the angle-integrated and differential cross 
section analysis. M^, cr^, /, a, and "radiative corr." describe the r] missing mass peak shape; D^f^g 
is the background missing mass function; other entries parameterize various detector properties. 
See text for details. 



with an energy-dependent width, Eq. ([8]) are used to fit the low W region. Various model 
calculations [321] in the past have found a small nonresonant contribution to the cross sec- 
tion and none is needed here. The simple shape describes the low W region well, but 
there are deviations for W > 1.6 GeV, presumably due to interference between S'ii(1535), 
511(1650), and nonresonant processes. Although the higher mass resonance is very near 
to the state we seek to describe, all analyses [9] find a very small rjN branching fraction 
for 5*11(1650). Therefore, we restrict the fit to W values less than 1.6 GeV. Two previous 
experiments 3, [l7| performed longitudinal/transverse separations in the late 1970's. Their 
results are consistent with no longitudinal component, albeit with large uncertainties. For 
the results presented here, the different beam energies have insufficient overlap in W and 
to separate these components. Under the assumption that the cross section is dominated 
by a single resonance and that Si is small, we can relate the Ai to the peak cross sections 

2 2 

extracted from the fit (see Eqs. (|1]),([5])): 

nl/2 



AliQ' 



R 



2mphri 



(12) 
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TABLE III: Binning details for the angle-integrated cross sections. For each bin, we show the 
minimum and maximum values of Q^, the energy of the electron beam for the data set, and the 
maximum value of W probed. The W bin width in all cases was 10 MeV. 
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1 SI 
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4.0 


2.30 


1.7 


2.1 


4.0 


2.30 


2.1 


2.5 


4.0 


2.13 


2.5 


2.9 


4.0 


1.93 


2.9 


3.3 


4.0 


1.72 



Consistent with Armstrong et al. [19| , a value of the full width of 150 MeV and an — > r]N 
branching ratio of 0.55 were used. The results of this determination of A\ are shown 

2 

in Fig. m along with some previous results converted to be consistent with our choice of 
Tr and 6^. The precise normalization of Ai depends on the choice of parameters for the 
contributing resonances, which are, as yet, not well determined. For instance, using the 
range of values listed in PDG for and 6^ leads to a 11% systematic uncertainty on Ai. 
While these uncertainties affect the absolute value of Ai, the shape of the dependence 

2 

is much better determined. More detailed understanding of this state is required to better 
determine absolute values of A\ and estimate the model dependence of those values. Using 
the choice of T and 6^, described above, our extracted values are consistent with the previous 
values at low Q'^ [3, 14, 16, [l7|, but with smaller uncertainties. At high Q^, there is 



moderate disagreement between the previously pubUshed results of Brasse et al. [18] (at 
= 2.0 and 3.0 GeV^) and Armstrong et al. [19] (at 2.4 and 3.6 GeV^). Our results match 
up nicely with Armstrong et al. and provide a precise determination of the shape of the 
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1.6 1.8 2.0 2.2 1.6 1.8 2.0 2.2 1.6 1.8 2.0 2.2 



W (GeV) 



FIG. 3: Angle-integrated cross sections {(y{'~lvP VP)) measured in all bins. Only statistical 
uncertainties are included in the data points. Systematic uncertainties are small compared to 
statistical errors and are not shown. The line represents the single Breit-Wigner fit. 

dependence of Ai from low up to their high determinations. 

The literature has various theoretical calculations of the photon coupling amplitude 
within the Constituent Quark Model (COM). Matching the slow falloff with has been 



difficult. We show two recent calculations 



40]. Aiello, Giannini, and Santopintc 



oQ 



use a 



hypercentral CQM and emphasize the importance of the 3-body quark force. Although this 
prediction gives the best agreement with our data of all the calculations, it falls off more 
rapidly with than the data. The Capstick and Keister calculation 2|] starts with the 
more traditional CQM, but use relativistic dynamics in a light-front framework. Although 
the two calculations use different approaches, the CQM is not well-defined and many other 
results are given in the literature. 
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FIG. 4: Extracted values for Ai for -yp 5ii(1535) transition for this experiment (filled squares). 

2 

There are overlapping data points at = 0.25 and 0.35 GeV^ coming from data taken with two 
different magnetic field settings. Filled circles show results from our previous smaller data set 
JLab results from Armstrong et al 



earlier publications 



11 



13 



14 



13, 



191 are shown as crosses. The open circles are results from 



181] . All results have been converted to represent a common 



width and branching ratio for S'ii(1535) ijp. The theoretical models are from constituent quark 
models of Capstick and Keister [2| (solid line) and Aiello, Gianinni, and Santopinto [40] (dotted 
curve) . 



B. Differential Cross Sections 



For larger bins in W and (see Table WV\ . we extract differential cross sections versus 
center-of-mass scattering angles of the rj {cos 6* and (p*). Each bin is labeled by its centroid. 
Results from this experiment are tabulated in the CLAS database [sol. For the = 0.8 
GeV^ bin, Figure [5] shows sample cross sections for four W bins. The first two W bins, 
W = 1.52 and 1.54 GeV are at the peak of the 5*11(1535) resonance. They show a dominant 
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TABLE IV: Binning details for the differential cross sections. For each bin, we show the 
minimum and maximum values of Q^, the energy of the electron beam for the data set, the 
maximum value of W probed, and the bin width in cos 0*. The W bin widths in all cases were 20 
MeV, while the bins were 30° wide. 
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2.5 


1.74 


0.2 
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2.1 


4.0 


2.00 


0.4 


2.1 


2.9 


4.0 


1.92 


0.4 



isotropic component due to the 5*11(1535) r]p process, but deviations from isotropy can 
be seen, especially at large 0*. BjW = 1.66 GeV, the non-isotropy is quite evident. The 
cross section falls monotonically as a function of cos^*, with the cross section for forward 
rj production consistent with zero. As W increases this feature changes dramatically. At 
W = 1.72 GeV, the forward-backward asymmetry of the distributions has reversed, with 
forward rj production favored, while backward production is close to zero. 

The ?7-MAID model j^], based on the MAID formalism |3o|, has been developed for r) 
electro- and photoproduction. This is an isobar model using a relativistic Breit-Wigner W 
dependence with form factors. Eight PDG 3* and 4* resonances of mass less than 1.8 GeV 
and nonresonant processes are included at the amplitude level. They fit the photoproduction 
data 



data 



12 



42l | and the dependence of the total cross section from electroproduction 
2l| in the 5*11(1535) region. The results of a calculation implementing this model 
are included in Figure [H These calculations roughly match the observed cross sections. 
However, the angular dependence predicted by r^-MAID does not agree with our data at W 
above the 5*ii(1535) region. The model was not fit to the differential cross sections of our 
previous work and the dependence of the higher mass resonances, e.g. Di5(1675), was 
taken from a quark model calculation rather than from data. 

For each W and bin, the differential cross sections are fit to a form that comes from an 
expansion of the response functions from Eq. ([2]) in terms of associated Legendre polynomials 
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FIG. 5: Sample differential cross sections for rjp in the center of mass frame for = 

0.8 GeV^ and selected W bins. Values for (p* symmetric about 180° have been averaged. (No 
information is lost this way; see Eq. ([2]).) Solid lines with an error band correspond to the resp onse 



function fit described in the text. Dashed lines correspond to the calculation of ry-MAID 
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(13) 



£=2 



The parameters Bg, Bi, and Ci depend on Q^, W, and e. They represent bilinear sums over 
contributing multipole amplitudes. Truncating to £ < 3, we determined the parameters in 
Eq. (fT3l) by a fit to the data. This truncation is motivated by three effects, i) The lightest 
known A^* resonance with £ > 3 is the Gi7(2190); if the dominant effects on the differential 
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cross section arise from interference with the dominant i = partial wave, terms above 
I = 3 should be negligible, ii) Fits to the r^-MAID predicted cross sections yield negligible 
contributions for terms higher than i = 3. iii) Good fits to the data are obtained with the 
truncated sum. 




W (GeV) 

FIG. 6: Results from fitting the angular distribution data of this experiment to Eq. (|13p . Coeffi- 
cients of the (f)* independent terms are shown, i.e. those that contribute to aT + eai- Contributions 
from both statistical and systematic sources are displayed. The dashed line is the r/-MAID predic- 
tion l29\ and the solid line is a four resonance fit to these terms. 

Results are shown versus W in Figures El El and [3 The quoted uncertainties contain 
both statistical and systematic uncertainties. We repeated the fit taking into account shifts 
in the cross section for each of the sources of systematic uncertainty studied in Section lllli 
The total systematic uncertainty on the extracted parameters is the sum in quadrature of all 
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FIG. 7: Same as Figured but showing the parameters corresponding to a lt- For the four resonance 
fit, these parameters are all zero. 

individual sources. We normalize our fitted Ai, A2, etc. to the isotropic term (Aq) in order 
to more clearly show the W and dependence of the shape of the differential cross section. 
For the ratios the resulting uncertainty is dominated by the uncertainty on the numerators. 

The isotropic component {Aq = atot/^'^) shows the same features as the angle-integrated 
cross sections: a dominant peak from the S'ii(1535) with additional structure above W = 1.6 
GeV. The other prominent term in the fit is Ai, which represents the slope of the differential 
cross section versus cos 6**. A structure in the W dependence of Ai was first seen in our 
previous publication 



2l| and is also seen by the GRAAL photoproduction experiment 42 1 



By examining the ratio Ai/Aq in the new data, we can study this structure in more detail. 
Two features stand out in this ratio: 

1. The ratio Ai/Aq is large and makes a rapid change from negative to positive values at 
W ^ 1.7 GeV. 



2. This structure is roughly independent of up to 2.5 GeV^. 
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FIG. 8: Same as Figure[6l but showing the parameters corresponding to utt- For the four resonance 
fit, these parameters are all zero. 

The simplest description for Ai is in terms of interference between S and P waves. In that 
case the rapid change in Ai between W = 1.66 and 1.72 GeV could be caused by one of the 
waves passing through a resonance. There are 2 P-wave resonances in this region, Pii(1710) 



and Pi3(1720). The former is rated 3* by the PDG 



to extract from data and is therefore controversial 43|. The latter is rated 4*, but is also 



, but its properties are very difficult 



poorly understood [lO|. Fits to CLAS tt+vt" electroproduction data [23] provided evidence 
that the existing baryon structure at W ~ 1.7 GeV should be changed. Their fits prefer 
either a greatly reduced pN decay branch for the existing Pi3(1720) resonance or a new I"*" 
state. In the present data, we cannot couple to a T = 3/2 state and are unable to distinguish 
between Pn and P13 states; we choose to use only a Pn state. If one describes the cross 
section using only Su and Pn partial waves, 

Ai _ 23?(E*+Mi_ ^ 

A' 



(14) 



|Eo+P + |Mi_|2- 

In this case, the rapid shift from backward to forward peaked cross sections would be due to 
a rapid change in the relative phase of the Eq+ and Mi_ multipoles because one of them is 
passing through resonance. The observation that this structure in Ai/Aq is approximately 
independent would then imply that Sn and Pn partial waves have a similar dependence. 
The values of shown in Fig. [7] are consistent with zero. These parameters measure 
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the aiT component of indicating that longitudinal amplitudes are not significant for this 
reaction (as was assumed in Section [IV Ap . The Ci parameters in Fig. [8] measure the cttt 
component of They are small indicating the components are also small for these 
values of Q^. 

To better understand the content of the ?7-MAID model, we also fit the parameters in 
Eq. f|T3|) to the predicted cross sections from that model. The extracted parameters are also 
included in Figures [H [TJ and [8] as dashed lines. The prediction has a broader Su peak than 
is seen in our data. Some structure is predicted in Ai arising from the Pii(1710), but the 
size of this effect is not nearly enough to match our data. At high W, ?7-MAID predicts a 
negative Ai in contrast to the significant positive value we observe. The model value of b^jN 
for Pii(1710) is much larger than the PDG value. Our data indicates the model value is 
incorrect. ?7-MAID contains many sources of D-wave contributions: Di3(1520), Di3(1700), 
and Di5(1675) in addition to nonresonant amplitudes. This produces a value for A2/A0 that 
matches our data for W < 1.6 GeV, where the Di3(1520) is the leading contribution. At 
larger W, agreement is poor. The model value of b^^N for -Di5(1675) is much larger than the 
PDG value; our data indicates this is incorrect. The prediction for A^/Aq is near zero, as 
are our measurements. 

Predictions for the air and utt terms are consistent with our measurements. C2 is the 
only term that is not negligible in //-MAID for our values of Q^. It arises from the A3/2 
amplitudes of the D-wave resonances interfering with the larger S-wave amplitude. Our 
data agrees with this general trend, but the effect is small compared to the uncertainties. 

To gain further understanding of the resonance content of our data, we did an additional 
fit to the differential cross section data using relativistic Breit-Wigner resonances according 
to Eqs. (!71)- (fT0l) . We fit the extracted parameters, up to W = 1.8 GeV, to a sum of four 
amplitudes for the following resonances: S'ii(1525), S'ii(1650), Pii(1710), and -Di3(1520). 
This set of resonances was determined empirically as the minimal set required to fit the 
general features of our data. Although the properties of Pii(1710) are very uncertain, it is 
an important contributor to this fit. We label it as Pn, but we cannot distinguish between 
Pii and Pi3(1720) in our data set; specifically, a P13 resonance would also give a rapid 
energy dependence in either A2 or C2 which we are unable to exclude with current statistical 
accuracy. Only the transverse response function, ctt is modeled, i.e. the Ai parameters. For 
the resonances with small contributions (5*11(1650) and Di3(1520)), we fixed the resonance 
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parameters to values obtained elsewhere. Masses and widths were set to average values 
from the Particle Data Group. For the Di3(1520) we used the dependence of r^-MAID. 
Following the assumption of the single quark transition model 4J|, the ratio of the strength 
of the 5*11(1650) to that of the 5*11(1535) was taken to be independent of . Motivated 
by the independence of A^j Aq in our data, we assumed the Pii(1710) had the same 
dependence as the 5*ii states. This left 12 variables in the fit: the masses and widths of 
the ^ii(1535) and Pii(1710), the relative strengths of the ^ii(1650) and Pii(1710) to that 
of the 5*11(1535), an overall strength of the Di3(1520), and the absolute strength of the 
5*11(1535) in each of the five bins. We view this as a simple fit. Our results should not be 
interpreted as a precise determination of resonance parameters, but rather as an indication 
of the dominant components needed in any future theoretical work. 

The results of this fit are also shown in Figure [61 The fit yields a reasonable, though 
not perfect description of our data. The isotropic term is described by the dominant 
5ii(1535) peak, modified by the smaller 5ii(1650). The deviation from a simple Breit- 
Wigner is described as a combination of destructive interference between the 5ii(1535) and 
5*ii(1650), and a small contribution from the Pii(1710). Including the 5*ii(1650) results in 
an extracted value of A\ for the 5ii(1535) which is 7% higher than that obtained with a 
single Breit-Wigner. The fitted width of the Pii(1710) is 100 MeV, which is consistent with 
the central (but very uncertain) PDG value. We cannot isolate the Pii(1710) photocoupling 
from that state's branching ratio into ?7p; we can only quote a ratio of ^ values (Eq. ([6])). 
The extracted value of ■Cirio/^isss is 0.22, which is about twice as large as in 77- MAID, and 
nearly an order of magnitude larger than that extracted from parameters of the Pii(1710) in 
the PDG. The -Di3(1520) primarily effects the quadratic term A2. Including this resonance 
is enough to give a reasonable description of the W dependence of A2. Our data do not 
require significant contributions from higher D-wsnq states present in the ?7-MAID model. 

We fit the structure in Ai/Aq with a smooth 5- wave and a rapidly changing P-wave. 
One could also describe this structure in terms of a new 5*-wave resonance interfering with 
P-wave component as in the model of Saghai and Li j^]. However, the amplitudes for the 
new resonance and the P-wave component must both fall off slowly with to reproduce 
the data. 

We also fit the (p dependence of the differential cross sections directly to Eq. ([2]) in order 
to obtain cr^ -|- ecr/,, ctlt, and axr as a function of W, and cos^^*. We choose to fit 
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the 0* dependence in terms of the parallel/perpendicular asymmetry (Asymrr) and the 
parallel/anti-parallel asymmetry (Asymi^)- 



where 



and 



Asym^nj. = 



+ (7_l' 
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FIG. 9: Extracted values for ax + eai, AsyrriLT, and Asym^^r as a function of cosO* for four 
selected W bins with = 0.8 GeV^. 
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For photoproduction, ai and ctlt don't contribute. A common polarization parameter 
is the parallel/perpendicular asymmetry, S. It is defined by Asym^-^ = eS; note also that 
S = ottI<^t- In electroproduction, the possible presence of a longitudinal term makes the 
relationships more complicated: 

Asymj,^, = ^^'^'^ (19) 

The data was analyzed in terms of these 3 response function combinations. A treatment 
of systematic uncertainties similar to that used for differential cross sections was applied. 
Figure [9] shows the values extracted from these fits for the same W and shown in Figure 
|5l Error bars display the systematic and statistical uncertainties. The quantity 0"^ + ecr^,, 
shows the same features we have discussed earlier. The values for Asym^r are consistent 
with zero in all distributions, but the size of the estimated error bars are a strong function of 
W . For W > 1.6, the total cross section is smaller than where the Su resonance dominates. 
Extraction of meaningful values for the dependence in this manner is therefore difficult. 



V. CONCLUSIONS 



Our extractions of Ai for the excitation of the 5*11(1535) cover a large range and match 
up well with Armstrong's results [19] at higher Q^. It should be noted again that there 
are significant model dependencies on describing the mass, width and branching ratio into 
rjp. These uncertainties lead to significant systematic uncertainties on the absolute scale of 
Ai. These uncertainties are common to all points currently determined, so the shape of the 
distribution is well determined. It becomes a significant challenge for theory to reproduce 
this shape. No existing model is able to describe the full range. 

Knowledge of the N* resonances in the region W ~ 1700 MeV is presently weak because 
the quality of older vrA^ tttiN and vrA^ — ^ rjN data is poor. The coupling of known 
P-wave resonances to rjN is thought to be very small. In this experiment, rapid energy 
dependence in the strength in the P-wave for coupling to rjN final states is found. With 
a simple resonance model, we are able to describe these data with significant coupling of a 
P-wave resonance to rjN. As with 5*11(1535), the falloff of this coupling must be very slow. 
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Although we can describe our measurements in terms of the Pii(1710), we cannot dis- 
tinguish between that or the Pi3(1720) with these data . Either resonance could produce 
the effect seen in Ai/Aq. The Pii(1710) is more poorly understood than the Pi3(1720), so 
it is easier to accommodate our data by altering the partial widths of the Pn rather than 
the P13. A large P13 could also produce effects in other terms. For instance, interference 
with a D-wave, would give a small contribution to As, but not significant compared to 

our uncertainties. A P13 resonance would have an A3 photo-excitation as well as Ai. Our 

2 2 

determination of Ai/Aq is sensitive to the Ai amplitude: a significant Aa amplitude could 

2 2 

also lead to large effects in cttt- 

Evidence for possible alterations in A^* P-wave resonances at masses about 1.7 GeV is 
accumulating. In addition to what is found in this experiment, double-pion production 
experiments in this same mass range 23|, |46|] are also unable to be described with models 



using existing information. Since different models are used to describe the different data 
sets, it is important to use a common model to describe the combined measurements from 
these (and other) reactions. Such a program may allow us to accurately determine the 
properties of the P-wave resonances in this region. 
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